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Abstract 

The aim of this paper is to solve the bispectral problem for bispectral oper- 
ators whose order is a prime number. More precisely we give a complete list of 
such bispectral operators. We use systematically the operator approach and in 
particular - Dixmier ideas on the first Weyl algebra. When the order is 2 the 
main theorem is exactly the result of Duistermaat-Griinbaum . On the other 
hand our proofs seem to be simpler. 

Introduction 

Bispectral operators have been introduced by F.A.Griinbaum (cf. [ pi| , p^ ]) in his 
studies on applications of spectral analysis to medical imaging. 

In the present paper we give complete classification of bispectral operators of 
prime order. We start with some definitions and results that are needed to state our 
results, as well as to make clear the connection with other research. 

An ordinary differential operator L{x,dx) is called bispectral if it has an eigen- 
function il){x,z), depending also on the spectral parameter z, which is at the same 
time an eigenfunction of another differential operator h.{z,dz) now in the spectral 
parameter z. In other words we look for operators L, A and a function ijj{x,z) 
satisfying equations of the form: 



f{z)^, 
e{x)i;. 



(0.1) 
(0.2) 



Although, as mentioned above, the study of bispectral operators has been stim- 
ulated by certain problems of computer tomography, later it turned out that they 
are connected to several actively developing areas of mathematics and physics - the 
KP-hierarchy, infinite-dimensional Lie algebras and their representations, particle 
systems, automorphisms of algebras of differential operators, non-commutative ge- 
ometry, etc. (see e.g. pHYl| , pHY^ , |BHY1 , Pw1 , |bWiI^W^ , Pg| ,|k|, |W1| , |^ ^ 
as well as the papers in the proceedings volume of the conference in Montreal ||BP| ). 

*E-mail; horozov@fmi.uni-sofia.bg and @math. bas.bg 
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In the fundamental paper |DG| Duistermaat and Griinbaum raised the problem 
to find all bispectral operators and completely solved it for operators L of order two. 
The complete list is as follows. If we present L as a Schrodinger operator 



L 



, d , 
' dx' 



+ u{x), 



the potentials u{x) of bispectral operators, apart from the obvious Airy {u{x) = ax) 
and Bessel {u{x) = cx~'^) ones, are organized into two families of potentials u{x), 
which can be obtained by finitely many "rational Darboux transformations" 

(1) from u{x) = 0, 

(2) from u{x) = — (j)^"^. 



Thus the classification scheme prompted by the paper |DG| is by the order of 
the operators. G. Wilson |W1| introduced another classification scheme - by the 
rank of the bispectral operator L (see the next section for definitions). In the 



above cited paper |Wl (see also Wilson gave a complete description of all 

bispectral operators of rank 1 (and any order). In the terminology of Darboux 
transformations (see [ BHY1[] ) all bispectral operators of rank 1 are those obtained 
by rational Darboux transformations on the operators with constant coefficients, i. 
e. L = p{dx) € C. In the above mentioned papers [DG, Wl] the classification is split 
into two, more or less independent parts. First, there is an explicit construction of 



families of bispectral operators of a given class (order 2 in | DG |; rank 1 in |W1[| ) 
The construction can be given in terms of Darboux transformations of "canonical" 
operators (a notion that needs clarification, see the last section for some comments). 
A second part should be to give a proof that, if an operator (in the corresponding 
class) is a bispectral one, then it belongs to the constructed families. 

In the last few years there has been increased activity | BIIY1| , |BHY5 , KRo, ^] 



in the the direction of constructing classes of bispectral operators (IpHYl, BIIY3 



KRc 



For a survey on this subject, see [BP, Hl| and the references therein. 



To the best of our knowledge, all known up to now families of bispectral operators 



can be constructed by the methods of | BHYl , BIIY3]. For a simplified exposition 
of these results, see the first part of | |H1| . A challenging problem is to prove that 
all the bispectral operators have already been found. A natural approach would 
be to divide the differential operators into suitable classes, e.g. - by order as in 



[ DG ] or by rank and to try to isolate the bispectral ones amongst them. In [ HM ] 
we have proposed another classification scheme - that is to consider the operators 
with a fixed type of singularity at infinity. The main result of that paper is the 
classification of bispectral operators possessing the simplest type of singularity at 
infinity - the Fuchsian one. My opinion is that all mentioned above classification 
schemes may help each other as seen from the main results here. 

In the present paper we return to the initial classification scheme - that of [DG[. 
We give a list of several families that contains all bispectral operators whose order 
is a prime number. Before stating the results we introduce some definitions and 
notations which will be used also throughout the paper. We are going to consider 
operators, normalized as follows: 
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N 

L = Y,yk{x)dt Vn = 1, FAr_i = 0. (0.3) 

fc=0 

It is well known that with the above normalization all the coefficients of L are 



rational functions (see | DG , Wl] or the next section). 

Now we can formulate the main result of the present paper. 

Theorem 0.1 An operator L, whose order is a prime number, is bispectral if and 
only if it belongs to one of the following sets: 

1. Generalized Airy operators: 

p-2 

A = dP + ajd^ - X, ttj G C; (0.4) 

2. Generalized Bessel operators: 

B = x-P{xd- (3i)...{xd- Pp), (3jeC; (0.5) 

3. Operators with constant coefficients: 

p-2 

C = dP + Y^ Ojd^ , Oj G C; (0.6) 
i=i 

4- Operators, obtained by monomial Darboux transformations from the Bessel 
operators having the property that at least one difference Pi — (3j G pZ, z ^ j; 

5. Polynomial Darboux transformations from operators with constant coeffi- 
cients. 



Remark 0.2 1. For p = 2 this is just the content of the classical result of Duistermaat- 



Griinbaum DG ]. Indeed, in that case we have Pi — P2 S 2Z. In DG ] this part 
of the theorem is formulated in a form close to this one. To obtain the potential 
u{x) = — (j)a;~^ mentioned above in 2) and corresponding to Pi = P2 = 1/2 we have 
to perform monomial Darboux transformation. Having in mind that a composition 
of monomial Darboux transformations is again monomial darboux transformation 
we get 2). The potential u{x) = from 1) corresponds to the only operator in 3., 
Theorem for p = 2. 

2. The different notions of Darboux transformations used here are explained in 
the next section. 



Theorem 3.1 is in fact a consequence of two slightly more general results. Their 



importance lies in the fact that they can be used for an induction process in further 



classification (see |II2|). In any case the proofs do not simplify when restricted to 
operators of prime order and seem to be more natural as performed here (see below). 

At the end we briefly review the organization of the paper. In Section 1. we recall 
some definitions and results, by now standard for the problem, with the purpose to 
fix the notation and the terminology. Section 2. treats the case of operators with 
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bounded (at infinity) coefficients Vj. We begin the section with an auxihary result. 
We show that from the normahzed operator L as in (|0.3| ) and satisfying only the 
so-called "ad-condition" with the polynomial 6 , the function f(z) can be chosen 
naturally and then one can build the operator A already normalized only in terms of 
L and 9. This will be needed after that to prove the following theorem for operators 
with bounded coefficients: 

Theorem 0.3 Let the rank of the bispectral operator L with bounded coefficients 
equals its order. Then it is a monomial Darboux transformation of a Bessel operator. 



The essential part of the proof is to establish the vanishing of the coefficients Vj at 
infinity. Then the result is contained in |HM]. 

In section 3. we consider operators with increasing coefficients. We first obtain 
a normal form for for their "leading terms" (subsection 3.1). Here we exploit once 
again (as in | HM[| ) crucial ideas of Dixmier analysis of the first Weyl algebra |Dx]. 
When the order is a prime number we get the non-vanishing (at infinity) part of L to 
be (generalized) Airy operator A. In the next subsection we develop some version of 
wave (pseudo-differential) operators, expanded in negative powers of Airy operators. 
This tool turns to be enough to prove (for operators of any order, not only prime) 
in subsection 3.3 the theorem: 



Theorem 0.4 Let the bispectral operator L = A +( vanishing at infinity perturba- 
tion). Then the perturbation is zero. 

We hope the careful reader has noticed some similarity with the analysis in |[DG[| for 
the latter class. We believe that our proof is simpler and more transparent (and for 
this reason works for higher order operators). One thing that we certainly benefited 



from [DG] is to realize that behind their calculations of the normal form of second 
order operators there stands Dixmier analysis on the first Weyl algebra Ai. On the 
other hand in both steps we use different approaches. 

The results of the present papers have been announced in the second part of the 



survey paper |H1[| . 
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of Bulgaria. 



1 Preliminaries 

In this section we have collected some terminology, notations and results relevant for 
the study of bispectral operators. Our main concern is to introduce unique notation 
which will be used throughout the paper and to make the paper self contained. 
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There are also few results which cannot be found formally elsewhere, but in fact are 
reformulations (in a form suitable for the present paper) of statements from other 
sources. 



1.1 

In this subsection we recall some definitions, facts and notation from Sato's theory 
of KP-hierarchy g |DJKM| , gwl needed in the paper. For a complete presentation 
of the theory we recommend also |Di, vM|. We start with the notion of the wave 
operator K{x,dx)- This is a pseudo-differential operator 

oo 

i^(x,9,.) = l + ^a,(x)9-^ (1.1) 
i=i 

with coefficients aj{x) which could be convergent or formal power (Laurent) series. 
In the present paper we will consider aj most often as formal Laurent series in x~^. 
The wave operator defines the (stationary) Baker- Akhiezer function ip{x,z): 

^{x,z)=Kix,d,)e'\ (1.2) 

Prom (|1.1| ) and (^]^) it follows that tp has the following asymptotic expansion: 

oo 

V'(x,z) = e^^(l + ^aj(x)z-J'), z^oo. (1.3) 
1 

Introduce also the pseudo-differential operator P: 

Pix,d,) = Kd,K-\ (1.4) 

The following spectral property of P, crucial in the theory of KP-hierarchy, is also 
very important for the bispectral problem: 

Ptp{x, z) = zip{x, z). (1-5) 

When it happens that some polynomial of P, say f{P), is a differential operator, we 
get that ip{x, z) is an eigenfunction of an ordinary differential operator L = f{P): 

Lij = f{z)i;. (1.6) 

It is possible to introduce the above objects in many different ways, starting with 
any of them (and with other, not introduced above). For us it would be important 
also to start with given differential operator L: 

L{x, d^) = + VN-2{x)d^-^ + ... + Vo{x). (1.7) 

One can define the wave operator K as: 

LK = Kf{d). (1.8) 
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An important notion, connected to an operator L is the algebra of operators 
commuting with L (see |Kr| , [BC|] ). This algebra is commutative one. The wave 
function ip{x,z) (defined in (|1.2|)) is a common wave function for all operators M 
from Al- 

Mi/j{x,z) =gM{z)ip{x,z). (1.9) 
We define also the algebra Al of all functions gM{z) for which ( |l.9| ) holds for some 



M ^ Al- Obviously the algebras Al and Al are isomorphic. Following [Ki] we 
introduce the rank of the algebra Al as the greatest common divisor of the orders 
of the operators m. Al- 

1.2 

Here we shall briefly recall the definition of Bessel wave function. Let /3 € be 
such that 

Eft^^^. (1.10) 

i=l 



Definition 1.1 [0, ||, |BHY1|] Bessel wave function is called the unique wave func- 



tion ^ p{x,z) depending only on xz and satisfying 

L^(x, dx)-^^{x, z) = z^-^p{x, z), (1.11) 
where the Bessel operator Lj3{x,dx) is given by (|0^). 
Because the Bessel wave function depends only on xz, ( |1.11| ) implies 

Dx^pix,z)=D,^pix,z), (1.12) 
Lp{z, d,)^pix, z) = x^^p{x, z). (1.13) 

Next we define monomial and polynomial Darboux transformations of Bessel 



operators. The definitions are slight modification of the definitions given in [ BHY1 |. 
Let h{Lp) be a polynomial in a Bessel operator. 

Definition 1.2 We say that the operator L is polynomial Darboux transformation 
of Lf} if there exist differential operators P{x, dx), Q{x, dx) and a polynomial h such 
that 

h{Lp) = QP (1.14) 
L = PQ (1.15) 

and the operator P{x,dx) has the form 

n 

P{x,dx) = x-^^Pk{x'')Dl (1.16) 

A:=0 

where pk are rational functions, p„ = 1. 
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We will use the following definition of monomial Darboux transformations. 



Definition 1.3 We say that the operator L is a monomial Darboux transformation 
of the Bessel operator iff it is a polynomial Darboux transformation with h{Lp) = 
L^, d G N. 



Remark 1.4 In |DG] the authors work with rational Darboux transformations. It is 
easy to show that a composition of rational Darboux transformations is a monomial 
Darboux transformation. 

We end this subsection by reformulating (in a weaker form) the main results, 
which we need from p3HYl , HM| ]. 



Theorem 1.5 The polynomial Darboux transformations of the Bessel operators are 
bispectral operators. 

Theorem 1.6 // a bispectral operator L has vanishing coefficients at infinity, it is 
a monomial Darboux transformation of a Bessel operator. 

1.3 



Here we recall several simple properties of bispectral operators following |DG, 
As we have already mentioned in the introduction we are going to study ordinary 
differential operators L of arbitrary order which are normalized as in ( p. 3]) , i.e. 
with V/v = 1 and Vat-i = 0. Assuming that L is bispectral means that we have also 
another operator A, a wave function iIj{x, z) and two other functions f{z) and 0{x), 
such that the equations ( |0.1| ) and ( |0^ ) hold. The following lemma, due to ]DG[ |, has 
been fundamental for all studies of bispectral operators. 

Lemma 1.7 There exists a number m G N, such that 

(adL)"+i0 = O. (1.17) 



For its simple proof, see [ DG| , |W1| . To the best of my knowledge this is the only 



property of bispectral operators that is used in their studies. It is widely believed 



that the condition ( 1.17 ), called the a.d- condition is equivalent to bispectrality (pro- 
vided (|0.3| ) holds). In what follows we assume only that the ad-conditon holds, i.e. 
we are not going to use the existence of f{z), A and ■0- For us it would be important 
to construct them only from L and 9 at least formally. This will be done in the next 
section. 

We will consider that m is the minimal number with this property. An important 
corollary of the above lemma is the following result. 



Lemma 1.8 Let the operator L be normalized as in ( |0.3| ). Then 
(i) The function 6{x) is a polynomial. 

(a) The coefficients aj in the expansion ( |1.1[ ) of the wave operator K are rational 
functions. 
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Proof. We repeat the simple proof following [ |W1U as the proof introduces some 
notions needed later. From the equation ( |1.17| ) it follows that 



{add^r+\K-^eK) = 0. 

On the other hand the kernel of the operator (ad 9^)™"*"^ consists of all pseudo- 
differential operators whose coefficients are polynomials in x of degree at most m. 
This gives that 

e{x)K = K&, (1.18) 
with a pseudo-differential operator : 

oo 

Q = Qo + Y,Qjdx^ (1-19) 

1 

whose coefficients Qj are polynomials of degree at most m. We have 6{x) = @q{x). 
This gives (i). Comparing the coefficients at dx'' we find that all the coefficients 
aj{x) of K are rational functions. □ 

Remark 1.9 We notice that at least one of the coefficients Qj has degree exactly 
m, where m from Lemma is minimal. This fact will be used later. 

The last lemma has as an obvious consequence one of the few general results, im- 
portant in all studies of bispectral operators. Noticing that the coefficients of L are 
polynomials in the derivatives of aj{x) we get 

Lemma 1.10 The coefficients of L are rational functions. 



Remark 1.11 Obviously the same results hold for the pair A, f{z), when imposing 
the conditions ( |0.3| ) on A. But here we need to derive this statement from conditions 
only on and L and a suitable choice of f{z) . This will be done in the next section. 

2 Operators with bounded coefficients 

In this section we are going to prove our main result for operators with coefficients 
bounded near infinity. As mentioned earlier we will consider slightly more general 
situation - operators for which the rank and the order coincide. The main part of 
the proof is to show that in this case the coefficients are in fact vanishing at infinity 



and hence the result follows from the main theorem in | HM | 
First I would like to fix the notation. Put 



p-2 

L = dP + ^Wj{x)d^, (2.1) 



j=0 



where Wj = cj + Vj{ constants and Vj{x) = 0{x ^). Next define the 

polynomial f{z) to be : 
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p-2 

f{z) = zP + Y,Cjz\ (2.2) 

j=0 

In this way we can rewrite our operator in the form 

p-2 

L = f{d) + Y,y,ix)d', (2.3) 

j=0 

where Vj{x) = 0{x~^). 

Our first goal is to to show that starting with the normahzed bispectral operator 
L the above choice of f{z) leads to a normalized operator A. The3 construction 
is formal, i.e. the wave function is a formal series, but this is enough for what 
follows. As all the auxiliary results (the three lemmas below) are slight modification 
of corresponding results from | HM[ | we omit their proofs. It is well known that one 



can present the operator L in the form: 

L = Kf{d)K-' (2.4) 



where the polynomial f{z) is defined in { \2.2\ j. In the next lemma, following [pG| we 
find the simplest restrictions on the coefficients of the wave operator K and on L. 



Lemma 2.1 (i) The coefficients Vj{x),j = N — 2, ...0 of L vanish at oo at least as 
(ii) The coefficients aj, j = 1, ... of the wave operator K vanish at least as 



X 



X 



This lemma allows us to introduce following | BHY4 | an anti- isomorphism h be- 
tween the algebra B of pseudo-differential operators P{x,dx) in the variable x and 
the same algebra B but in the variable z. More precisely B consists of those pseudo- 
differential operators 



p = X^p,(x-l)5-^ 

k 

for which there is a number n G Z (depending on P) such that all expressions 
x"'pj{x~^), j = k, k + 1, . . . are formal power series in x~^. The involution 

b:B — >B' 

is defined by 

oo 

6(p)e^^ = Pe^^ = z-^pj{d~^)e''^, for PeB (2.5) 

k 

i. e. b is just a continuation of the standard anti-isomorphism between two copies 
of the Weyl algebra. In what follows we will use also the anti-isomorphism 

bi:B — > B', bi{P) = 6(AdxP). (2.6) 
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Obviously b and 61 can be considered as involutions of B and without any ambiguity 
we can denote the inverse isomorphisms b~^,b^^ : B' — > B by the same letters. 
Since the operators K and = K~^6K are from B we can define two operators S 
and A as follows: 



S{z,d,) = b{K{x,d,)), (2.7) 
A{z,d,) = 6(G). (2.8) 



Explicitly one has 



and also 



S = Y^ z-^a,{d,) = J2 ajiz)d-\ ao = 1 (2.9) 
i=o j=o 

00 m 

A{z,d,) = Y^z-^e,{d,) = Y,Hz)dl (2.10) 

j=0 j=0 



where 7^ (see Remark 1.9) and the coefficients Aj and aj should be viewed as 
formal power series. As in |HM] we can prove 

Lemma 2.2 The coefficients aj of the operator S are rational functions. 

From the last lemma it follows that A is normalized as required in (^]^). Denote 
temporarily by r the degree of the polynomial 6, i. e. if 0{x) = + . . .. 

Lemma 2.3 With the choice of f(z) as in ( p.2| ) the coefficients Aj of the operator A 
are rational functions and A satisfies ( p.2D . The degree of 9 r = m and 



Am = 1, Am-1 = 0. (2.11) 

The point in the last lemma is that the normalization of A is a consequence 
of the suitable choice of the polynomial f{z). Of course the wave function is only 
formal but this suffices for the proof of the main theorem. 

Now we are ready to give the classification of operators with bounded near 
infinity coefficients of operators with the same rank and order. 

Let us fix the notation. Put 

N-2 

L = d^ +^Vj{x)d^, (2.12) 

j=0 

where Vj = cj + Wj{x), cj are constants and VF(x) = 0{x~^). As before define the 
polynomial f{z) to be : 

N-2 

f{z) = z^+^c,z^. (2.13) 

j=0 

With this choice of f{z) as we know (see ) the operator A is with rational coefficients. 
Now we are ready to prove the main result of this section. 
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Theorem 2.4 // the rank of the operator L with bounded at infinity coefficients 
equals its order N then all constants Cj = 0, i.e. the coefficients Vj = 0{x~^). 

Proof. We are going to use again the ad-condition ( |1.3| ). Choose the number m so 
that a(i^^^(A) 7^ and ad'J'^^{A) = 0. Simple computation shows that 

adf^^^iA) = i-irmlif {z)r (2.14) 
On the other hand we have 

ad^(e) = Q^0 (2.15) 

Obviously 

[L,Q]=0. (2.16) 

This gives that Q E A^. From this and from the fact that the rank of L is equal to 
its order N we get that Q is a polynomial in L: 

Q = qrV + qr-iU-^ + . . . , G C, (2.17) 

where the coefficient qr 7^ 0. Using the involution hi we get 

hi{Q) = qrnz) + Y,<ljf'i^) (2-18) 



Using (2.14) we have the following string of identities: 

61 (Q) = 6i(aC(0)) = (-l)™(ad7(,)(A)) = m!(/(z))-. (2.19) 
In this way (|1|) and (p^ yield 



r-l 

Irfiz) + Y.qjfHz) = m\{f'{z)r. (2.20) 

j=o 

Now we are going to compare both the degrees and the first two coefficients of 
the two hand-sides of ( 2.20| ). First notice that comparing the degrees of the leading 
terms gives: 

{N - l)m = rN. 

This gives that m is divisible by A^, i.e. m = sN and r = s{N — 1). Next, comparing 
the coefficients at the highest degree we get qr = m\N^\ Suppose that some of the 
coefficients Cj of f{z) are not zero . Denote the second non-zero term after by 
CkZ^. Computing the coefficient at the second non-zero terms in both sides ( p.2[)| ) 
we get k = N — 1. This is a contradiction to the normalizing condition cn-i = 0. □ 
From the above theorem and from the main result in [ HM| (see also Section 1.2, 



Theorem nM) we get the proof of Theorem 0.3. 
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3 Operators with increasing coefficients 



3.1 Normal forms 

Let L be a bispectral operator, normalized as in (|0.3| ), i.e. 

N~2 

L = d^+Y,yji^)di- (3.1) 

j=0 

We will consider that for some j the corresponding coefficient Vj is increasing at 
infinity, i.e. has Laurent expansion at infinity of the form: 

Vj{x)= J2 (3.2) 



m=—oo 



where rj > and aj^rj 7^ 0. We call rj the order of Vj. In what follows we are 
going to use several properties shared by the first Weyl algebra Ai and the larger 



algebra TZ[d] of differential operators with rational coefficients. Following [Dx| we 
define ffitration in Tl[d]. Let p, o" G M. We put wt{x) = p, wt{d) = a. Next the 
weight of L is given by the following definition: 

Definition 3.1 Assume that L = K5" + K-i5""-^H Wo is an arbitrary element 

of TZ[d]. For each term V{x)dl. define its weight 

Vp,„{V{x)di) = p{oTdV) + ai. 

Then the number 



Vp^a{L) := max Vp^a{Vi{x)d'') 

0<i<n 

will be called {p, c7)-order of L. 

The second definition associates to each differential operator from TZ[d] a (p, a)- 
homogeneous polynomial. 

Definition 3.2 Assume the notation of the previous definition and denote by I{L) 
the set {i S {0, 1, • • • , n}|fp^o-(^i'9*) = Vp^a-iL)}- The polynomial / G C[x,x^^,y] 
defined as: 

/(x,y) = ^a,x-'i^>y^ (3.3) 

where Oj G C are uniquely determined from the expansion 

Vi = Ojx"'^'^^' + {lower order terms), 

will be called polynomial associated with L. The operator Lq =: f{x,dx) '■ will 
be called the homogeneous part of L. Here, as usual the columns :: denote normal 
ordering, i.e. the differentiation is pushed to the right. 
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Consider again the bispectral operator (3.1). In what fohows we assume that p 



and a are positive integers. It is always possible to choose them in such a way that 
the polynomial / associated with L has at least two terms of the kind 

f = yN + ax^yP + ..., (3.4) 

where A; > 0, p > 0. For this purpose one can use the Newton polygon. More 
precisely denote by E{L) the set of points {m,j), such that am,j 7^ 0, where am,j 
is from ( |3.2| ). Consider the plane with the points of E{L) and take the convex 
closure of E{L) (the Newton polygon). Then draw a line passing through the point 
{0,N) £ E{L) and another point, say (fc,j) G E{L) with A; > and such that 
the Newton polygon remains below the line. This line is unique - all other points 
(A;i,ji) € E{L) with the same property lie on it. Then one can find a non-zero 
solution in integers of the equation Na = kp + ja. In our situation both p and a 
are positive as j < N and at least one Vj is icreasing. Notice that with this choice 
the polynomial f{x,y) £ C[x,y]. 

Denote also by g{x) = x^ the polynomial associated to 9{x). Our goal is to find 
severe restrictions on the polynomial /. We are going to use the following result 
which is a slight modification of a particular case of the fundamental Proposition 



7.3 from [Dx|. 



Lemma 3.3 Suppose the element L G TZld] acts on a element G nilpotently, i.e. 

ad2'(G) =0. m > 1 

Let f and g be the polynomials associated to L and G and f contains at least two 
terms. If p and a are positive integers and Vp^a{L) > p + a, then one of the following 
cases holds: 
(a) 

r = 9', (3.5) 

where s and r are the weights of g and f , 
(h) a > p, p divides a and 

f = X"(X" + fiY)\ (3.6) 

(c) p > a, a divides p and 

f = Y'^iY'^ + ^X)^ (3.7) 

(d) p = a and 

f = {Y + \xY{Y + ^iXf (3.8) 
where n>0, A;>0, m>0 and ^, A € C. 

Remarks on the proof of the lemma. The proof essentially repeats that of lemma 
7.3 from |Dx|. We mention the minor differences. While in [Dx] all the polynomials 



(p, (T)-associated with the elements of the Weyl algebra belong to C[x,y] , here we 
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work in C[x,x ^,y]- The fact that the polynomials belong to C[x,2/] is needed in 



[ Dx | mainly to speak about their roots in some algebraic closure of C(y) (respec- 
tively - C(x)),when considered as polynomials in x (respectively in y). But the ring 
considered here has the same property. 

Lemma 3.4 Suppose that k > 1 (k is from ( p.4| ) ). Then either Vp^cr{L) > p + a or 
L cannot act nilpotently on 6. 

Proof. Writing ( p.4| ) in the form / = y^{y^~^ + ax^) + . . . we notice that one can 
take p = N — p and a = k > 2 First suppose N = k = 2 and p = 0. Then L 
cannot be nilpotent (see |Dx| ). We give a slightly different proof here. One can 
easily see that if Vp^^iG) = 21 > 0, then Vp^aadLiO) = 21 unless f^ = g^. But this 
is not possible as the polynomial associated with adL(^) contains a term with y of 
degree one. Denote by gm the polynomial associated to ad™(0). By induction on m 
we see that the resulting applications of adi always contain terms where the power 
of y is less than 2. This shows that the equation f^ = g"^ is impossible. 

Now we can suppose that either p > or max(A^, k) > 2> (we recall that both 
> 2 and k > 2). In the first case we have Vp^^iL) = Na>kN>N + k> 
N -p + k> p + a. Ifp>0 then Vp^a{L) = Na = Nk>N + k>p + a. The strict 
inequality is a consequence of the fact that max{N, k) > 3. □ 

Lemma 3.5 Let L act on nilpotently. The polynomial f has the form 

f = {y^-x)\ r>l. (3.9) 

Proof. If the term with highest power k in (|3.4| ) is 1 then / has the form / = 
y^{y^ — Ax). In the case when k > 2 from Lemma f3.3| we know that / = y^{y^ — Xx)^ 
or / = (y — Xx)°'{y — fix)^, a + (3 = k, i.e. we have one of the cases b), c) or 
d) as = g^ is impossible. By applying the automorphism (^)-i (here 

= ^) ^r,fi{y) = y + A'^'") to / and to we reduce the last case to the 
previous ones, i.e. we assume that / = y^{y^ — Xx)^, where n > 1 and k > 1. 
Without loss of generality we can assume that A = 1. Our goal will be to show 
that if n > 1 and k > 1 then ad£ cannot be zero for any s G N. We will show 
that at least some of the terms with highest weight will be preserved. Suppose that 
n > 1, A; > 1. First we continue the automorphisms from Ai to its skew- 
field. Notice that they preserve the filtration. For this reason we are going to work 
only with the homogeneous part as the the rest of the terms have no impact on it. 
Apply the automorphism ^r,i of to / and to x'. Using that $r,i(9) = d and 
^r,iix) = x + d'^ we obtain 



^rAd'^id' -x)^) = (-1)'=9V, 



$^,i(x') = {x + d^'Y. 

(x + a^)' = ^4x^9"('-^) + ..., 



Consider now adfg„^fc-,(x + 9^)'. Write 
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where by ... we denote the lower weight terms. Then by hnearity 

I 

adlQ„^,^{x + = c'jadlQ^^,^{x^d'^'-^^) + ... (3.10) 

j=0 

Let k > n. Simple computation gives that 

a4„,.(x') = [Hinl + {k- n)i]9^("-i)]x'+^('=-^) + . . . 
j=0 

Asn>l, />1, k — n > the coefficient at the term of highest power in x 
is positive for any s > 1, which shows that ( |3.1[1| ) cannot be zero for any s. Now 
suppose that n > k > 1. Consider 

s-l 

adl^^,id'n = [{-ly l[[lrk + {n- k)jW^''-''>+'^ x^''-^^ + ... 

j=0 

By the same argument the coefficient at the highest power in d is not zero for any 
s. This shows that either n = or A; = 0. But from the assumption (|2.7D it follows 
that k cannot be zero. □ 

Remark 3.6 Note that the above result gives a normal form for the leading terms 
of all bispectral operators with increasing coefficients of any order. 

Now assume that the order of L is a prime number. This gives that k = 



1, N = r. Using Lemma 3.5 we obtain that: 
Lemma 3.7 The operator L has the form 

N-2 N-2 



L = d^ +Y, ajd^ -^+Y1 Wj{x)d^, (3.11) 

j=l j=0 



where lim^.^oo Wj{x) = and aj € C. 
We will call the operator 

N-2 



the principal part of L. Following the terminology of |BHY1| A is the (generalized) 
Airy operator. 



3.2 Airy PDO's 
Let 



N-2 

A = d^ ajd^ - X (3.12) 
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be the generalized Airy operator. Our aim here is to develop a calculus of pseudo- 
differential operators written in terms of inverse powers of Airy operators in complete 
analogy with the standard one, described in sect. 1.1. All the results of the section 
are obtained for any order A^, i.e. without assuming that the order is prime. Let 
^{x) be a nonzero function in KerA, i.e. 

A^{x) = 0. (3.13) 

Then the function ^(x, z) = (f){x + z) satisfies the equations: 



A{x, dx)'i'{x, z) = z^{x, z), (3.14) 
A{z, d^)^{x, z) = x^{x, z). (3.15) 

Obviously the function ^'(x,z) satisfies also the equation 

d^^{x,z)=d,^{x,z). (3.16) 

The equations ( |3.14| )-( |3.16| ) define an anti-involution b on the Weyl algebra Ai (see 
[ pHY2| , [WIU ), acting on the generators A,dx of Ai by 



biA{x,d,)) = z, (3.17) 

b{d,) = d,. (3.18) 

From (KM one easily finds that 

b{x) = A{z,d,), (3.19) 

which would be used later. 

Now define the algebra Bi of pseudo-differential operators of the type: 

oo 

P{x,d^)= aj{x,d,)A-^, (3.20) 

j=-m 

with operator coefficients of the form: 

N-l 



= ^3,k{x)d\ (3.21) 



fc=0 

where the functions aj^kix) are formal Laurent series: 

oo 

a,-fc(x)= (3.22) 

s=r 

and the index r depends only on P (but not on j!). Then the anti-automorphism b 
can be continued on Bi as 
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oo N—1 oo oo 

b{P)= z-^ Y^P^p'^A-^ = Ybsiz,dM-'{ 

j=-m k=0 s=r s—r 

Introduce the "wave operator" K as follows: 

oo 

K = 1 + Ymj{x,dr,)A-\ 



z . 



(3.23) 



where 



N-l 



k=0 

Let L be a differential operator of the form: 

L = A^ + Vi^iA^'^ + ..., 



where 



N-l 



(3.24) 
(3.25) 

(3.26) 
(3.27) 

t 

(3.28) 

The coefficients aj^kix) from ( 3.26| ) of the expansion of K can be found by induction 
from the equation: 

LK = KA^ (3.29) 

Multiplying the above equation from the right by A,A^^... one computes the coef- 
ficients aj^k{x) in the expansion ( 3.24| ) - (3.26) of K in terms of the functions Vj^^. 
We would particularly be interested in the case / = 1. 

In what follows we assume that the operator L is bispectral, it satisfies equation 
of the form (p. 11), together with an equation of the form (|0.2|). In that case as in 



fc=0 



Then one can find a wave operator (not unique) K of the form ( |3.24| ) so that 



L = KA^K-^. 



[ DG , Wl | one can prove the following lemma. 

Lemma 3.8 The coefficients aj^k of the operator K are rational functions. 



Proof. We mimic the well known proof (see [DG, |W1| ). Write 

{^dLTiO) = 0. 

{ad{A^)T{K-^eK). 



This is equivalent to 
Put 



Q = K-^OK = YejA-\ 

j=0 



(3.30) 
(3.31) 

(3.32) 



where 
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7V-1 



k=0 

This gives 

iadA^nOj) = 0. 
The leading terms of the above equation give 



'j,N-l'~' -\- ■ ■)'^ -t- . . . -h 



(0]7 . . + . . . = 

Here in the brackets containing the coefficients at s = 2, . . . , N the 

dots after %,Ar-s denote expressions of derivatives of 6j,M-r-, r < s oi order not lower 
than m. Then obviously by induction we get that ah (9j^Li = 0, n = 1, . . . , A^, 
which shows that they are polynomials of degree d < m. Then the computation of 
the coefficients aj{x) of K is performed as in [DG, [Wl| (see also sect. 2). We see 
that they are rational functions. □ 

3.3 Proof of the main theorem for operators with increasing coef- 
ficients 

Let L be an operator of order A'^ with Airy principal part 
A = d^ + aN-2d^-^ + ... + aid-x, 

i.e. 

Af-2 
j=0 

and 

lim Vj{x) = 0. 

X— >oo 



Using the techniques of the previous subsection we present L in the form 

L = A + V = KAR-^ (3.33) 

Our goal will be to show that bispectrality, and in particular - the rationality 
of the coefficients aj^k{x) in the expansion (|]2|) and (|3^ ) of K, implies that the 



perturbation V{x,dx) = 0, which is equivalent to aj^k = for all j,k. But first we 
need some notation and auxiliary results. 
From the equation 
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LK = KA, (3.34) 

we can compute recursively the coefficients rrij of the operator K. For this we will 
need some formulas to compare the coefficients of the two sides of ( |3.34| ). Introduce 
the operators bj, Cj, Uj, Wj, j = 1, 2, . . . by: 



[A,mj] = bjA + Cj] (3.35) 
V{x, d)mj = UjA + Wj. (3.36) 

We will need to order the monomials in m and related to it expressions as follows: 

Definition 3.9 We say that the monomial x^^d^^ is of higher order than the mono- 
mial x^^d^^ if ri > r2 or ri = r2 and ki > k2- We will call the number r the height 
of m, if the highest order term of m is of the type c.x^d^, c ^ 0. We denote this 
number by ht{m). 

In other words we use lexicographic ordering in the set of the monomials x^d^ but 
the order is only the power of x. 

In what follows we are going to use the abbreviation l.o.t. (for lower order terms) 
compared to some operator m with the meaning that they are lower than at least 
one of the terms in m. We will need also the following lemma: 

Lemma 3.10 In the above formulas we have 
(i) 

N-l 

bj = Na'j^kd'"'^ + l-o-t. (3.37) 

k=l 

N-l 

Cj = Naj Qd^-^ + ^ {xNaj^j, + kaj^k)d^-^ + Let. (3.38) 
k=i 

(a) ht{cj) = ht{bj) + 1. 

Proof. The proof is straightforward computation. To avoid two indices we will 
suppress the dependence on j (it is irrelevant at that moment) . We have 

N~l N~l 



Ao m = m{ akd^) — xm + 

Nat,d^+^-^ + l.o.t. (3.39) 

k=l k=0 
N~l N-l 

mo A = <^kd^) - xm + - ^ kakd^~^ . (3.40) 



k=l k=0 

Subtracting (|3lO| ) from ( |3^ we get 



N-l N-l 



[A, m] = Y, Na'^d^-^^'^ + ^ "^kd^^^ + l.o.t.. (3.41) 



fc=0 fc=0 
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Split the first sum into two parts as follows. One of them contains derivatives from 
to 2N — 2; the second will contain the rest of the them. Then we have for the 
first part 



N-l 



N-1 



k=0 k=l 

Next use the identity = A — J2k=i o-k + x to get 

N-l 
k=l 

N-l N-l 

Naod^-^ + ( ^ Na[d''-^)A + x( ^ Na'^d''-^) + ... 



k=l 



k=l 



where the dots represent terms with derivatives of a^. Then we repeat the same 
procedure to all terms (including the ones in l.o.t. from (|3.41| )) containing with 
k > N. After finite number of steps we get ( 3.37 ) and (|3.38| ). The second part of 
the lemma follows immediately from the first one. □ 



Lemma 3.11 (i) The equation ( |3.34| ) is equivalent to the equat 



ions: 



bi + V + Ui = 

bj+i + Cj + Wj + Uj+i = 0, J = 1, . . . 
(a) The coefficients ofV and bi behave at infinity as x'"^ 



(3.42) 
(3.43) 



Proof. The first part is simply comparing the coefficients. Indeed, writing in detail 
(|]3|) we get 

A + AmiA~^ + ... + V + VmiA~^ + ... = A + mi + ... 
Using (l3!35|) and we can simplify the last equation to 

V + bi + Ui + ... = 0, 

where . . . denote the purely pseudo-differential part. This gives dll). Multiplying 
( |3.34| ) by A, A"^, etc. from the right and arguing in the same manner we get ( p.43|) . 

To prove the second part we use (|3!3^ ) with j = 1 and (|1|). Notice that the 
leading terms of bi are derivatives of a rational functions. Being equal to the leading 
terms of V they vanish. Hence they vanish at least of order x~'^. □ 



Lemma 3.12 The following inequalities hold: 

ht{Uj) < ht{bj) 
ht{Wj < ht{cj) 



1, 
1. 



(3.44) 
(3.45) 
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Proof. The proof is similar to that of Lemma 



3.10 



Using ( p. 36 ) we obtain 



2N-1 N-l 



2N-1 



fc=l s=k 



k=0 



For k = N, . . . , N — 1 put Vj^k = Uj^k-N- As above spht the sum into two parts, the 
first one containing the terms with d^, k > N: 



N-l 



N-l 



Vmj = Uj^kd')d^ + E V,,kd' 



k=0 



k=0 



Again use the identity 9^ = A — ^^^i ak + x several times to get the first sum in 
the form: 



7V-1 



N-l 



k=0 

Then obviously we have: 



k=0 



N-l 



Uj = Uj^kd^ + l.o.t. 

k=0 

N-l N-l 



(3.46) 
(3.47) 



fc=0 



fc=0 



Now using that the order at infinity of F is x ^ we get that ht{Uj) < ht{m) — 2, 
ht{Wj) < ht{mj) — 1. Prom the last inequalities we get ( 3.44 ) and ( 3.45 ). □ 



Now we are ready to finish the proof of the Theorem 3.4. 
Proof, of Theorem We recall that we have to show that V = 0. Assume that 
some of the coefficients Vj are not zero. Then we shall compute the leading terms of 
the operators bj recursively using ( |3.43| ) and Lemma p. 11 and taking into account 



the estimates ( 3.44| ) and ( 3.45| ) First notice that that the highest order term in bi is 
of the type aix'^^d'', ai 7^ 0, with k < N and si < 0. Suppose that the highest order 
term in bj is ajx'^^d'', aj 7^ 0. Then the highest order term in bj^i is computed, using 



(Q), ( |337D and ( p8| ) to be aj+ix'^^+^d'' with aj+i = -aj{N{sj+l)+k)/N{sj+l). 



Having in mind k < N we get that a^+i 7^ 0. After a finite number of steps we will 
get that for some j the corresponding Sj = —1. But this contradicts the fact that 
the highest order term of bj is a derivative of a rational function. □ 



4 Final remarks on the proof and comments 
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4.1 Proof of Theorem 0.1 



Essentially we already have performed the proof of the main theorem. We just have 
to notice that when the order of L is prime and there are coefficients increasing 



at infinity Lemma |3.7| and Theorem |0.4| give that the operator is Airy and hence 
bispectral. 

If the coefficients of L are bounded (at infinity) and the order is prime, then 
using the fact the the rank divides the order we get that either the rank of L is 1 



or it is equal to its order. The latter case is treated in Theorem |0.3| . If the rank is 



one then this is the main result of |W1|. Finally the inverse part, i.e. that all the 



operators listed in Theorem 3.1 are bispectral is the main result of |BHY1|.Q 



4.2 Comments 

Here I would like to make some speculations on eventual continuation of the classi- 
fication. It seems to me that the methods of [lBHY3|| (see also [HI] for more details) 
will be enough to construct all bispectral operators. Assuming that then the clas- 
sification should be: 1) find all "basic" bispectral operators; 2)show that Darboux 
transforms reduce any bispectral operator to a "basic" one. 

Having in mind the constructions in |DG, BHYl, KRo, |W1| it seems natural to 



consider basic those operators that have as few singularities as possible and generate 
their centralizers. Then in view of the main result of [ [H2| one class of operators that 
certainly should be considered as "basic" is the class of bispectral operators L in 
the Weyl algebra that together with some other operator Q satisfy the "canonical 
commutation relation" (CCR): 



[L,Q] 



1 



(4.1) 

all bispectral operators in 



More precisely they are basic because according to 
the Weyl algebra are simply polynomials in operators L that satisfy (4J.). It is 
tempting to believe that the CCR is enouhg for bispectrality: 



Conjecture 4.1 If the operators L and Q satisfy the CCR (4T) then they are 
bispectral. 



This conjecture seems to be a difficult one as it is easily shown (cf. 
equivalent to the famous conjecture of Dixmier-Kirillov: 



1]) to be 



Conjecture 4.2 If the operators L,Q satisfy the CCR (4.1) then they generate 
the Weyl algebra Ai. In other words, any endomorphism of the Weyl algebra is an 
automorphism. 



The Bessel operators as well as other examples from |BHY3| show that one needs 
also to add to the list of the basic operators L the ones that together with some 
other operator Q satisfy the following "string" equation 

[L,Q]=L (4.2) 

Our, maybe insufficient experience, suggests that the above two classes contain all 
basic bispectral operators. A more precise conjecture is the following one: 
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Conjecture 4.3 Every bispectral operator is a polynomial Darboux transformation 
of a bispectral operator satisfying either ( [4.1D or (|4.2| ). 



To make the classification of the basic bispectral operators more explicit intro- 
duce the following notation. Denote by the algebra spanned by a Bessel operator 
Lfy, and D. Then 



Conjecture 4.4 (1) A bispectral operator satisfying (4.1) belongs the Weyl algebra. 



(2) A bispectral operator satisfying (4.2) belongs to one of the algebras Ba- 



Some progress could be achieved if one finds analogs for B^ of the results from 
[H2]. One difficulty would be to extend Dixmier results for these algebras. 

All the above conjectures aim to a complete classification of bispectral operators. 
A more modest goal is to find the bispectral operators with some properties. Here 
is a conjecture in this direction. 

Conjecture 4.5 An operator L with bounded at infinity coefficients is bispectral if 
and only if it is a polynomial Darboux transformation of a Bessel operator. 



This conjecture seems natural in view of the "if part, obtained in |BHY1|. It will 
be very useful either to give a proof of Wilson's result about rank one operators 
without using algebraic-geometric arguments or modify his proof in higher rank 
situation. Despite of the many interesting results for higher rank solutions of KP- 
hierarchy I have not found a construction suitable for our purposes. In any case the 
above conjecture seems to me to be within the reach of the existing tools unlike the 
previous ones. 
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